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This  study  investigates  the  long-run  and  the  causal  relationship  between  greenhouse  gas  emissions, 
economic  growth  per  unit  of  energy  use  and  energy  consumption  in  Pakistan  by  using  cointegration 
theory  and  Granger  causality  test.  In  addition,  this  study  compares  the  influencing  magnitude  of 
greenhouse  gas  emissions  on  energy  consumption  by  employing  variance  decomposition  techniques 
over  a  36-year  time  period,  i.e.,  between  1975  and  2011.  The  study  employed  the  Johansen  cointegration 
technique  to  estimate  the  long-run  relationship  between  the  variables,  Granger  causality  test  to  find  the 
influencing  directions  while  variance  decomposition  was  used  to  compare  the  influencing  magnitude 
between  the  variables.  The  study  was  limited  to  a  few  variables,  including  greenhouse  gas  emissions 
(such  as  agricultural  methane  emissions;  agricultural  nitrous  oxide  emissions;  C02  emissions  and 
combustible  renewables  and  waste);  GDP  per  unit  of  energy  use  and  energy  consumption,  in  order  to 
manage  robust  data  analysis.  Finding  suggests  that  energy  consumption  acts  as  an  important  driver  for 
increase  in  greenhouse  gas  emissions  in  Pakistan.  The  results  indicate  that  on  average,  causality  runs 
from  energy  consumption  to  greenhouse  gas  emissions  but  not  vice  versa.  Energy  consumption  does  not 
Granger  cause  agricultural  nitrous  oxide  emissions  via  both  routes  which  confirm  ‘neutrality  causal 
hypotheses'  between  the  variables.  Greenhouse  gas  emissions  are  closely  associated  with  economic 
growth  per  unit  of  energy  use  and  energy  consumption  in  Pakistan.  Variance  decomposition  analysis 
shows  that  among  all  the  greenhouse  gas  emissions,  combustible  renewables  and  waste  exerted  the 
largest  contribution  to  changes  in  energy  consumption  in  Pakistan.  The  results  suggest  that  consuming 
energy  by  industries  is  becoming  more  increasing  by  every  coming  day,  but  it  signifies  the  fact  that 
emissions  caused  by  consuming  the  energy  is  causing  harm  to  the  society  by  enlarge,  and  it  reaches  to 
the  conclusion  where  we  are  heading  toward  a  systematic  fatal  disorder  of  our  society. 
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1.  Introduction 

Many  chemical  compounds  found  in  the  Earth’s  atmosphere  act 
as  "greenhouse  gases.”  These  gases  allow  sunlight  to  enter  the 
atmosphere  freely.  When  sunlight  strikes  the  Earth's  surface,  some 
of  it  is  reflected  back  towards  space  as  infrared  radiation  (heat). 
Greenhouse  gases  absorb  this  infrared  radiation  and  trap  the  heat 
in  the  atmosphere.  Over  time,  the  amount  of  energy  sent  from  the 
sun  to  the  Earth's  surface  should  be  about  the  same  as  the  amount 
of  energy  radiated  back  into  space,  leaving  the  temperature  of  the 
Earth's  surface  roughly  constant.  Many  gases  exhibit  these  “green¬ 
house”  properties.  Some  of  them  occur  in  nature  (water  vapor, 
carbon  dioxide,  methane,  and  nitrous  oxide),  while  others  are 
exclusively  human-made  like  gases  used  for  aerosols  [1], 

One  way  of  attributing  greenhouse  gas  (GHG)  emissions  is  to 
measure  the  embedded  emissions  of  goods  that  are  being  con¬ 
sumed.  Emissions  are  usually  measured  according  to  production, 
rather  than  consumption  [2].  At  the  global  scale,  the  key  green¬ 
house  gases  emitted  by  human  activities  are: 

•  Carbon  dioxide  (C02)  -  fossil  fuel  use  is  the  primary  source  of 
C02.  The  way  in  which  people  use  land  is  also  an  important 
source  of  C02,  especially  when  it  involves  deforestation.  Land 
can  also  remove  C02  from  the  atmosphere  through  reforesta¬ 
tion,  improvement  of  soils,  and  other  activities. 

•  Methane  (CH4)  -  agricultural  activities,  waste  management,  and 
energy  use  all  contribute  to  CH4  emissions. 

•  Nitrous  oxide  (N20)  -  agricultural  activities,  such  as  fertilizer 
use,  are  the  primary  source  of  N20  emissions. 

•  Fluorinated  gases  (F-gases)  -  industrial  processes,  refrigeration, 
and  the  use  of  a  variety  of  consumer  products  contribute  to 
emissions  of  F-gases,  which  include  hydrofluorocarbons  (HFCs), 
perfluorocarbons  (PFCs),  and  sulfur  hexafluoride  (IPCC,  [36]). 


Global  greenhouse  gas  emissions  further  break  down  by  the 
economic  activities  that  lead  to  their  production  (see,  Fig.  1 ). 

Ensuring  that  natural  resources  are  consumed  and  waste  is 
produced  at  sustainable  rates  represent  major  contemporary 
challenges  for  the  world.  Recognition  of  these  challenges  resulted 
in  the  endorsement  in  2000  of  environmental  sustainability  as  one 
of  the  Millennium  Development  Goals  (MDGs)  to  be  achieved  by 
2015.  However,  by  2003  global  rates  of  consumption  and  waste 
production  were  estimated  to  be  at  least  25%  higher  than  the 
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Fig.  1.  Global  greenhouse  gas  emissions  by  source. 

Source:  IPCC  [36], 


capacity  of  the  planet  to  provide  resources  and  absorb  waste  and 
this  rate  may  have  risen  as  high  as  50%  by  2007  [3]. 

Climate  change  is  an  area  that  has  become  increasingly 
important  in  recent  years  and  raises  issues  of  global  justice  and 
equity,  whereas  the  richer  industrialized  countries  are  primarily 
responsible  for  greenhouse  gas  emissions;  it  is  the  poorer  devel¬ 
oping  countries  that  would  most  heavily  bear  the  costs  of  climate 
change.  Global  climate  change  leads  to  an  increased  number  of 
weather-related  disasters  such  as  floods  and  droughts,  which 
cause  food  shortages  and  famine.  However,  agriculture  not  only 
suffers  from  environmental  problems,  but  it  also  contributes  to 
them,  through  pollution,  overgrazing,  and  release  of  greenhouse 
gases  [4]. 

The  quality  of  the  natural  environment  is  not  only  an  extremely 
important  issue  from  the  point  of  view  of  individual  survival  but  it 
will  also  emerge  as  one  of  the  principal  human  security  issues  in 
Pakistan  [5],  Pakistan's  agrarian  economy  is  heavily  dependent  on 
river  water  provided  by  melting  glaciers.  Pakistani  cities  are  facing 
problems  of  urban  congestion,  deteriorating  air  and  water  quality 
and  waste  management  while  the  rural  areas  are  witnessing  rapid 
deforestation,  biodiversity  and  habitat  loss,  crop  failure,  desertifi¬ 
cation  and  land  degradation.  The  major  climate  change  threats  to 
Pakistan  are; 

•  Considerable  increase  in  frequency  and  intensity  of  extreme 
weather  events,  coupled  with  erratic  monsoon  rains  causing 
frequent  and  intense  floods  and  droughts. 

•  Increased  siltation  of  major  dams  caused  by  more  frequent  and 
intense  floods. 

•  Increased  temperature  resulting  in  enhanced  heat-  and  water- 
stressed  conditions,  particularly  in  arid  and  semi-arid  regions, 
leading  to  reduced  agricultural  productivity. 

•  Further  decrease  in  the  already  scanty  forest  cover  from  too 
rapid  change  in  climatic  conditions  to  allow  natural  migration 
of  adversely  affected  plant  species. 

•  Increased  health  risks  and  climate  change  induced  migration,  etc. 

The  above  threats  are  the  cause  of  major  survival  concerns  for 
Pakistan,  particularly  in  terms  of  the  country's  water,  food,  and 
energy  security  considerations  [6], 

There  is  widespread  agreement  among  energy  economists  and 
policy  analysts  that  the  ever-increasing  energy  demand  is  the 
major  contributing  factor  to  anthropogenic  greenhouse  gas  emis¬ 
sions.  Energy  is  essential  to  all  economic  activities,  however,  the 
increasing  attention  given  to  global  warming  and  climate  change 
has  renewed  spur  to  find  the  relationship  between  environmental 
pollutants,  energy  consumption,  and  economic  growth.  Nowadays, 
mitigation  assessment  of  greenhouse  gas  emissions  has  become  an 
integral  part  of  the  national  and  international  climate  policy 
agenda  [7]. 

The  rise  in  global  energy  demand  has  raised  questions  regard¬ 
ing  energy  security  and  increased  the  focus  on  diversification, 
generation  and  efficient  allocation.  The  answer  lies  in  the  attain¬ 
ment  of  optimal  energy  mix  through  fuel  substitution  by  promot¬ 
ing  energy  efficiency  and  renewable  energy  and  interregional  co¬ 
operation  [6],  Pakistan's  economy  has  been  growing  at  an  average 
growth  rate  of  almost  3%  for  the  last  four  years  and  demand  of 
energy  both  at  production  and  consumer  end  is  increasing  rapidly. 
Table  1  shows  the  supply  of  energy  in  Pakistan  over  a  period  from 
2001-2002  to  2010-2011. 
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Table  1 

Primary  energy  supply  and  per  capita  availability. 
Source :  GoP  [6], 


Year  Energy  supply  Per  capita 


Million  TOE 

Change  % 

Availability  (TOE) 

Change  % 

2001-2002 

45.07 

1.5 

0.32 

-1.25 

2002-2003 

47.06 

4.4 

0.31 

0.00 

2003-2004 

50.85 

8.1 

0.34 

6.25 

2004-2005 

55.58 

9.3 

0.36 

5.88 

2005-2006 

58.06 

4.5 

0.37 

2.78 

2006-2007 

60.62 

4.4 

0.38 

2.70 

2007-2008 

62.92 

3.8 

0.39 

2.63 

2008-2009 

62.55 

-0.6 

0.38 

-2.56 

2009-2010 

63.09 

0.9 

0.36 

-5.6 

2010-2011 

64.52 

2.3 

0.36 

0.00 

Energy  needs  are  indelibly  linked  to  Pakistan's  economic  and 
sustainable  growth  capabilities.  Pakistan  has  been  in  increasing 
demand  across  the  various  areas  of  energy  sources.  With  a 
growing  economy  and  the  desire  for  vast  production  and  con¬ 
sumption  across  the  country,  the  energy  demands  remain  high. 

The  objective  of  this  study  is  to  empirically  investigate  the 
impact  of  energy  consumption  on  greenhouse  gas  emissions  in 
Pakistan.  The  more  specific  objectives  are: 

I.  To  estimate  whether  there  is  a  long-run  relationship  among 
greenhouse  gas  emissions,  GDP  per  unit  use  of  energy  and 
energy  consumption  in  Pakistan. 

II.  To  explore  the  influencing  directions  between  greenhouse  gas 
emissions  and  energy  consumption. 

III.  To  compare  the  influencing  magnitude  of  greenhouse  gases  on 
energy  consumption  in  Pakistan. 

The  study  is  arranged  in  the  following  manner.  An  introduction 
has  been  discussed  in  this  section.  Section  2  shows  related 
literature  review.  Data  source  and  methodological  framework  are 
mentioned  in  Section  3.  Results  are  shown  in  Section  4.  Final 
section  concludes  the  study. 


2.  Literature  review 

There  is  now  a  wide  and  growing  body  of  literature  on  the 
potential  impacts  of  global  warming.  According  to  [8], 

“Scientific  research  on  greenhouse  impacts  has  so  far  almost 
entirely  concentrated  on  the  benchmark  case  of  warming  under 
an  atmospheric  C02  concentration  of  twice  the  pre-industrial 
level  (2  x  C02).  As  a  consequence  studies  on  the  economic  costs 
of  global  warming  have  tended  to  concentrate  on  the  same 
benchmark”. 

Lee  [9]  explores  whether  energy  conservation  policies  can  be 
implemented  in  countries  with  the  same  level  of  development,  for 
this  purpose,  author  employs  Granger  non-causality  testing  pro¬ 
cedure  to  re-investigate  the  relationship  between  energy  con¬ 
sumption  and  income  in  11  major  industrialized  countries  over  a 
period  of  1960-2001.  The  results  do  not  support  the  view  that 
energy  consumption  and  income  are  neutral  with  respect  to  each 
other,  except  in  the  case  of  the  United  Kingdom,  Germany  and 
Sweden  where  a  neutral  relationship  is  found.  Bi-directional 
causality  in  the  United  States  and  unidirectional  running  from 
energy  consumption  to  GDP  in  Canada,  Belgium,  the  Netherlands 
and  Switzerland  are  found.  According  to  [10], 


“Generally,  gas  emissions  due  to  materials  or  product  transpor¬ 
tation  and  production  are  considered  part  of  local  energy 
consumption  and  are  attributed  to  the  system  or  country 
where  the  fuels  are  distributed  without  considering  why  these 
fuels  are  burnt,  who  benefits  from  their  use  and  thus  who  is 
responsible  for  their  consumption”. 

Zhang  and  Cheng  [11  ]  investigate  the  existence  and  direction  of 
Granger  causality  between  economic  growth,  energy  consump¬ 
tion,  and  carbon  emissions  in  China  over  the  period  of  1960-2007. 
The  results  suggest  a  unidirectional  Granger  causality  running 
from  GDP  to  energy  consumption,  and  a  unidirectional  Granger 
causality  running  from  energy  consumption  to  carbon  emissions 
in  the  long  run.  Soytas  and  Sari  [12]  investigate  the  long  run 
relationship  between  economic  growth,  carbon  dioxide  emissions 
and  energy  consumption  in  Turkey  over  a  period  of  1960-2001. 
The  results  show  that  carbon  emissions  seem  to  Granger  cause 
energy  consumption,  but  the  reverse  is  not  true.  Soytas  et  al.  [13] 
investigate  the  effect  of  energy  consumption  and  output  on  carbon 
emissions  in  the  United  States  over  a  period  of  1960-2004.  The 
results  reveal  that  income  does  not  Granger  cause  carbon  emis¬ 
sions  in  the  US  in  the  long  run,  but  energy  use  does. 

Shahbaz  et  al.  [14]  examine  the  linkages  among  economic 
growth,  energy  consumption,  financial  development,  trade  open¬ 
ness  and  C02  emissions  over  the  period  of  1975Q1-2011Q4  in  the 
case  of  Indonesia.  The  empirical  findings  indicate  that  economic 
growth  and  energy  consumption  increase  C02  emissions,  while 
financial  development  and  trade  openness  compact  it.  The  caus¬ 
ality  analysis  has  shown  the  feedback  hypothesis  between  energy 
consumption  and  C02  emissions.  Economic  growth  and  C02 
emissions  are  interrelated,  i.e.,  bi-directional  causality  while, 
financial  development  Granger  causes  C02  emissions.  Marcotullio 
et  al.  [15]  identify  and  explore  greenhouse  gas  emissions  from 
urban  areas  in  Asia  at  the  regional  level;  and  further  explore 
covariates  of  urban  greenhouse  gas  emissions.  The  results  suggest 
that  urban  areas  account  between  30%  and  38%  of  total  anthro¬ 
pogenic  greenhouse  gas  emissions  for  the  region  and  that  emis¬ 
sion  per  capita  averages  from  urban  areas  are  lower  than  those  at 
the  national  level.  Important  covariates  for  total  urban  greenhouse 
gas  emissions  include  population  size,  density  and  growth  rate, 
income  per  capita,  development  status  and  elevation. 

Al-mulali  et  al.  [16]  investigated  the  long  run  relationship 
between  urbanization,  energy  consumption  and  carbon  dioxide 
emission  in  seven  regions,  namely,  East  Asia  and  Pacific,  East 
Europe  and  Central  Asia,  Latin  America  and  the  Caribbean,  Middle 
East  and  North  Africa,  South  Asia,  Sub-Saharan  Africa,  and  Western 
Europe  over  a  period  of  1980-2008.  The  results  show  that  while 
84%  of  the  countries  have  a  positive  long  run  relationship  between 
urbanization,  energy  consumption,  and  carbon  dioxide  emission, 
only  16%  the  countries  have  mixed  results.  Some  countries  have  a 
negative  long  run  relationship  and  others,  especially  low  income 
countries  have  no  relationship  between  urbanization,  energy 
consumption,  and  carbon  dioxide  emission.  Furthermore,  a  one 
way  long  run  relationship  between  energy  consumption  and 
carbon  dioxide  emission  and  urbanization  was  found  in  a  number 
of  countries  while  a  one  way  long  run  relationship  between 
urbanization  and  energy  consumption  and  carbon  dioxide  emis¬ 
sion  was  found  in  other  countries.  Zaman  et  al.  [17]  investigate  the 
influence  of  agricultural  technologies  on  carbon  emissions  in 
Pakistan  by  using  annual  data  from  1975  to  2010.  The  results 
reveal  that  agricultural  technologies  act  as  an  important  driver 
for  increase  in  carbon  emissions  in  Pakistan.  Results  indicate 
that  unidirectional  causality  runs  from  agriculture  machinery  to 
carbon  emissions  but  not  vice  versa.  Agricultural  technologies  are 
closely  associated  with  economic  growth  and  carbon  emissions  in 
Pakistan. 
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Shahbaz  et  al.  [18]  investigate  the  relationship  between  energy 
(renewable  and  nonrenewable)  consumption  and  economic 
growth  in  case  of  Pakistan  over  the  period  of  1972-2011.  The 
results  confirm  cointegration  between  renewable  energy  con¬ 
sumption,  nonrenewable  energy  consumption,  economic  growth, 
capital  and  labor  in  case  of  Pakistan.  The  findings  show  that  both 
renewable  and  nonrenewable  energy  consumption  add  in  eco¬ 
nomic  growth.  Capital  and  labor  are  also  important  determinants 
of  economic  growth.  The  VECM  Granger  causality  analysis  vali¬ 
dates  the  existence  of  feedback  hypotheses  between  renewable 
energy  consumption  and  economic  growth,  nonrenewable  energy 
consumption  and  economic  growth,  economic  growth  and  capital. 
According  to  Ali  and  Nitivattananon  [19,  p.  77] 

“Pakistan's  industrial  and  residential  sectors  are  vibrant  con¬ 
sumers  of  energy  and  C02  emitters  among  all  other  sectors  of 
the  city.  Sparse  trees  in  the  city  and  reduced  agriculture  areas 
by  more  than  one-half  in  2009  compared  with  those  in  1975 
are  the  main  reasons  for  increased  energy  use  and  reduced  C02 
emissions  from  agriculture  sector  as  well.  However,  all  the 
other  sectors  have  increased  their  C02  emissions  in  an  escalat¬ 
ing  trend”. 

Ali  et  al.  [20]  assess  the  expected  future  effects  of  national  and 
global  biofuels  policies  on  agriculture  markets  and  food  prices  in 
Pakistan.  The  results  show  that  by  2020,  global  mandates  on 
biofuels  will  significantly  affect  the  prices,  production  and  trade 
of  major  feedstock  crops  such  as  sugarcane,  maize,  soybean  and 
rapeseed,  especially  in  the  USA,  Brazil  and  EU.  Global  biofuels 
developments  are  projected  to  increase  the  prices  of  maize,  rape- 
seed,  soybean,  and  sugarcane  in  Pakistan.  According  to  Zaman 
et  al.  [21,  p.  13] 

“Environmental  impacts  are  the  unwanted  byproduct  of  eco¬ 
nomic  activities.  Inadvertently,  humans  alter  environmental 
conditions  such  as  the  acidity  of  soils,  the  nutrient  content  of 
surface  water,  the  radiation  balance  of  the  atmosphere  and  the 
concentrations  of  trace  materials  in  food  chains”. 

Sustainable  development  has  become  a  major  policy  goal  over 
recent  decades,  but  this  goal  is  still  very  often  seen  as  conflicting 
with  traditional  economic  and  energy  policy  goals.  More  recently, 
the  cost  of  climate  change  has  received  more  attention  and  an 
economic  case  for  sustainable  development  has  been  made. 
Investigations  on  the  links  between  energy  consumption  and 
sustainable  development,  by  contrast,  are  still  in  their  infancy, 
except  perhaps  for  the  potential  of  ‘green’  jobs  [22].  Recognizing 
the  knowledge  gap  on  the  linkages  between  energy  consumption 
and  sustainable  development,  this  study  initiate  to  examine  the 
impact  of  energy  consumption  on  greenhouse  gas  emissions  in 
Pakistan  and  discusses  the  implications  for  policy  makers  and  puts 
forward  recommendations  to  be  considered  by  the  policy 
community. 


3.  Data  source  and  methodological  framework 

The  annual  time  series  data  is  employed  for  the  Pakistan 
economy  over  a  period  of  1975-2011.  All  relevant  data  is  taken 
from  World  Bank  development  indicators  which  are  published  by 
World  Bank  [23].  Green  house  gas  emissions  (GHG)  (such  as 
agricultural  methane  emissions;  agricultural  nitrous  oxide  emis¬ 
sions;  C02  emissions  and  combustible  renewables  and  waste)  and 
electric  power  consumption  are  used  to  shed  light  on  the  possible 
impact  of  energy  consumption  on  green  house  gas  emissions  in 
Pakistan  in  the  presence  of  GDP  per  unit  of  energy  use.  Table  2 
shows  the  list  of  the  variables  and  their  expected  causalities 
between  them. 
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ENRG 


Fig.  2.  Data  trend  for  selected  variables. 

Source :  World  Bank  [23], 

All  these  variables  are  expressed  in  natural  logarithm  and 
hence  their  first  differences  approximate  their  growth  rates.  The 
data  trends  are  available  for  ready  reference  in  Fig.  2. 

Fig.  3  highlights  in  schematic  fashion  the  methodological  approach 
adopted  in  this  study.  According  to  this  framework,  green  house  gas 
emissions  have  been  checked  on  energy  consumption  through  GDP 
per  unit  of  energy  use. 

Following  four  equations  (Panel  A-D)  are  used  to  assess  the 
impact  of  green  house  gases  on  energy  consumption  in  Pakistan,  i.e„ 


Panel  A  :  ENRG,  AGRMET,  GDPENRG  ( 1 ) 

Panel  B  :  ENRG,  AGRN1T,  GDPENRG  (2) 

Panel  C  :  ENRG,  C02,  GDPENRG  (3) 

Panel  D  :  ENRG,  CRWASTE,  GDPENRG  (4) 


where  ENRG  represents  energy  use  (kt  of  oil  equivalent);  AGRMET 
represents  agricultural  methane  emissions  (thousand  metric  tons 
of  C02  equivalent);  AGRNIT  represents  agricultural  nitrous  oxide 
emissions  (thousand  metric  tons  of  C02  equivalent);  C02  represents 
carbon  dioxide  emissions  (metric  tons  per  capita);  CRWASTE  repre¬ 
sents  combustible  renewables  and  waste  (metric  tons  of  oil  equiva¬ 
lent)  and  GDPENRG  represents  GDP  per  unit  of  energy  use  (constant 
2005  PPP  $  per  kg  of  oil  equivalent). 


CRWASTE 


Y  ears 


3.1.  Econometric  framework  of  the  study 

The  time  series  data  often  show  the  property  of  non- 
stationarity  in  levels  and  the  resulted  estimates  usually  provide 
spurious  results.  Thus,  the  first  step  in  any  time  series  empirical 
analysis  was  to  test  for  presence  of  unit  roots  to  remove  the 
problem  of  inaccurate  estimates.  The  other  important  step  was  to 
check  the  order  of  integration  of  each  variable  in  a  data  series  in 
the  model  to  establish  whether  the  data  under  hand  suffer  unit 
root  and  how  many  times  it  needed  to  be  differenced  to  gain 
stationarity. 

The  test  for  cointegration  consists  of  two  steps:  first,  the 
individual  series  are  tested  for  a  common  order  of  integration.  If 
the  series  are  integrated  and  are  of  the  same  order,  it  implies 
cointegration.2  Dickey  and  Fuller  [24,25]  devised  a  procedure  to 
formally  test  for  non-stationarity.  The  Augmented  Dickey  Fuller 
(ADF)  test  is  used  to  test  the  stationarity  of  the  series.  The  ADF  test 
is  a  standard  unit  root  test:  it  analyzes  the  order  of  integration  of 
the  data  series.  These  statistics  are  calculated  with  a  constant,  and 
a  constant  plus  time  trend,  and  these  tests  have  a  null  hypothesis 
of  non-stationarity  against  an  alternative  of  stationarity. 


2  If  the  series  are  integrated  with  the  mixture  of  order  of  integration  i.e.,  I 
(0)  and  7(1 ),  it  implies  bonds  testing  approach  which  was  proposed  by  Pesaran  et  al. 
137]. 
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Fig.  3.  Research  framework. 
Source:  self  extract. 


Johansen’  cointegration  tests  applied  on  the  series  of  same 
order  of  integration,  i.e„  /  (1)  series  which  determine  the  long  run 
relationship  between  the  variables.  When  series  are  co-integrated 
of  order  1,  trace  test  (Johansen's  Approach)  indicates  a  unique  co¬ 
integrating  vector  of  order  1  and  hence  indicates  the  long  run 
relationship,  in  the  multivariate  case,  if  the  /(I )  variables  are  linked 
by  more  than  one  co-integrating  vector,  the  Engle-Granger  [26] 
procedure  is  not  applicable.  The  test  for  cointegration  used  here  is 
the  likelihood  ratio  put  forward  by  Johansen  and  Juselius  [27], 
indicating  that  the  maximum  likelihood  method  is  more  appro¬ 
priate  in  a  multivariate  system.  Therefore,  this  study  has  used  this 
method  to  identify  the  number  of  co-integrated  vectors  in  the 
model.  The  Johansen  and  Juselius  method  has  been  developed  in 
part  by  the  literature  available  in  the  field  and  reduced  rank 
regression,  and  the  co-integrating  vector  Y  is  defined  by  Johansen 
as  the  maximum  Eigen-value  and  trace  test  or  static,  there  is  Y  or 
more  co-integrating  vectors.  Johansen  [28]  and  Johansen  and 
Juselius  [27]  proposed  that  the  multivariate  cointegration  metho¬ 
dology  could  be  defined  as 

(ENRG)t  =  (GHG,  GDPENRG) 

where  GHG  represents  green  house  gas  emissions,  which  is  a 
vector  of  elements.  Considering  the  following  autoregressive 
representation: 

K 

ENRGt  =  n°  +  2  ffj(ENkC)t-i  + 

T  =  l 

Johansen's  method  involves  the  estimation  of  the  above  equa¬ 
tion  by  the  maximum  likelihood  technique,  and  testing  the 
hypothesis  H0;  of  “r”  co-integrating  relationships,  where 

r  is  the  rank  or  the  matrix  n(OzrzP),  Y  is  the  matrix  of  weights 


with  which  the  variable  enter  co-integrating  relationships  and  £  is 
the  matrix  of  co-integrating  vectors.  The  null  hypothesis  of  non¬ 
cointegration  among  variables  is  rejected  when  the  estimated 
likelihood  test  statistic  ej>A  =  -nZf=r+1ln(l-2;)}  exceeds  its  critical 

value.  Given  estimates  of  the  eigen-value  (2,)the  eigen-vector  (£,) 
and  the  weights  (¥',),  we  can  find  out  whether  or  not  the  variables 
in  the  vector  (ENRG)  are  co-integrated  in  one  or  more  long-run 
relationships  among  (GHG,  GDPENRG). 

This  study  investigates  the  influence  of  Pakistan's  greenhouse 
gases  on  energy  consumption  from  two  perspectives.  One  is  to 
conduct  the  modified  Granger  causality  and  Johansen  cointegra¬ 
tion  tests  to  explore  the  influencing  directions  between  green¬ 
house  gases  and  energy  consumption,  respectively;  the  other  is  to 
compare  the  influencing  magnitude  of  greenhouse  gases  on 
energy  consumption,  based  on  the  vector  error  correction  model 
(VECM)  and  variance  decomposition  approach. 

In  order  to  undertake  the  modified  version  of  Granger  causality 
for  a  VAR  model  with  3  lags  (k= 2  and  dmax=l).  we  estimate  the 
following  system  of  equations: 


"ENRG 

"ENRG 

"ENRG 

AGRMET 

=  Aq  +  Ai 

AGRNIT 

+  a2 

C02 

GDPENRG 

GDPENRG 

GDPENRG 

"ENRG 

"fit" 

+^3 

CRWASTE 

+ 

£2 1 

GDPENRG 

£3  c_ 

where  A i  and  A3  are  the  3x3  matrices  of  coefficients  with  A0 
being  a  3  x  1  identity  matrix,  and  et  are  the  disturbance  terms  with 
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zero  mean  and  constant  variance.  From  Eq.  (5)  we  can  test  the 
hypothesis  that  Pakistan's  greenhouse  gases  does  not  Granger 
cause  energy  consumption  with  the  following  hypothesis,  i.e., 


where  a\2  are  the  coefficients  of  the  greenhouse  gases  in  the  first 
equation  of  the  system  presented  in  Eq.  (5).  Besides,  we  can  test 
the  opposite  causality  from  Pakistan's  energy  consumption  to 
greenhouse  gases  in  the  following  hypothesis: 


where  are  the  coefficients  of  the  energy  consumption  variable  in 
the  second  equation  of  the  system  presented  in  Eq.  (5).  It  should  be 
noted  that  we  incorporate  the  variable  GDP  per  unit  of  energy  use 
in  to  Eq.  (5)  to  avoid  the  omitted  variable  bias  when  we  examine 


the  Granger  causality  bias  when  we  examine  the  Granger  causality 
between  greenhouse  gas  emissions  and  energy  consumption. 


4.  Results  and  discussions 

4.3.  Cointegration  among  greenhouse  gases  and  energy  consumption 

The  present  study  performs  the  augmented  Dickey-Fuller 
(ADF)  unit  root  tests  for  variables  with  regard  to  their  stationary 
properties.  The  detailed  results  are  shown  in  Table  3. 

The  results  reveal  that  all  variables  in  this  study  are  non¬ 
stationary  at  their  level,  though,  stationary  at  their  first  differ¬ 
ences.  We  say  that  all  variables  are  integrated  of  order  one,  i.e.,  I 
(1)  series  at  1%  level.  Fig.  4  shows  the  data  trend  at  their  first 
difference. 


Table  3 

Augmented  Dickey-Fuller  (ADF)  test  on  the  levels  and  on  the  first  difference  of  the  variables  (1975-2011). 


Variables  Level 


First  difference 


Decision 


Constant 

Constant  and  trend 

Constant 

Constant  and  trend 

ENRG 

0.954  (0) 

-3.035  (0) 

-5.449*(0) 

-5.437*(0) 

Non-stationary  at  level  but  stationary  at  first  difference,  i.e.,  /  (1) 

AGRMET 

-2.144  (1) 

-0.579  (0) 

-4.450*(0) 

-5.045*(0) 

Non-stationary  at  level  but  stationary  at  first  difference,  i.e.,  /  (1) 

AGRNIT 

1.294  (0) 

-2.903  (0) 

-6.336*(4) 

-6.265*(4) 

Non-stationary  at  level  but  stationary  at  first  difference,  i.e.,  /  (1) 

co2 

-0.122  (0) 

-3.012  (1) 

-6.577*(0) 

-6.438*(0) 

Non-stationary  at  level  but  stationary  at  first  difference,  i.e.,  /  (1) 

CRWASTE 

0.375  (0) 

-2.167  (0) 

-4.992*(0) 

-5.049*(0) 

Non-stationary  at  level  but  stationary  at  first  difference,  i.e.,  /  (1) 

GDPENRG 

-2.412  (0) 

-3.233  (0) 

-4.438*(0) 

-4.359*(0) 

Non-stationary  at  level  but  stationary  at  first  difference,  i.e.,  /  (1) 

Note:  The  null  hypothesis  is  that  the  series  is  non-stationary,  or  contains  a  unit  root.  The  rejection  of  the  null  hypothesis  is  based  on  MacKinnon  [39]  critical  values,  i.e.,  at 
constant:  -3.639,  -2.951  and  -2.614  are  significant  at  1%,  5%  and  10%  level,  respectively.  While  at  constant  and  trend:  -4.252,  -3.548  and  -3.207  are  significant  at  1%,  5%  and 
10%  level,  respectively.  The  lag  length  is  selected  based  on  SIC  criteria,  this  ranges  from  lag  zero  to  lag  four. 


D  (AGRMET)  D  (AGRNIT) 


D  (ENRG) 


D  (GDPENRG) 


Fig.  4.  Data  trend  at  first  difference  for  selected  variables. 
Source:  author's  own  efforts  with  the  help  of  EViews  version  6. 
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After  that,  we  take  energy  consumption  (ENRG)  as  the  depen¬ 
dent  variable  and  each  greenhouse  gases  and  GDP  per  unit  use  of 
energy  (GDPENRG)  together  as  the  independent  variables,  and 
then  the  Johansen  cointegration  among  them  is  tested  according 
to  Johansen  [28],  From  the  results  in  Table  4,  we  find  that  except 
C02,  all  other  greenhouse  gas  emission  have  at  least  one  coin¬ 
tegration  relationship  with  energy  consumption  at  5%  level. 
Therefore,  we  may  say  that,  greenhouse  gases  have  significant 
long-term  equilibrium  with  energy  consumption  in  Pakistan. 


4.2.  Causality  among  greenhouse  gases  and  energy  consumption 

Subsequently,  we  conduct  the  modified  Granger  causality  tests 
by  Toda  and  Yamamoto  [38]  for  greenhouse  gas  emissions  and 
energy  consumption.  The  variable  GDP  per  unit  use  of  energy 
(GDPENRG)  is  incorporated  as  an  explanatory  variable  to  avoid  the 
omitted  variable  bias.  Results  are  shown  in  Table  5. 

The  results  reveal  that  “energy  consumption  (ENRG)  does  not 
Granger  cause  agriculture  methane  gas  emissions  (AGRMET)’’  is 
rejected  at  5%  level  but  not  vice  versa.  This  result  suggests  that 
there  is  unidirectional  causality  runs  from  energy  consumption  to 
agricultural  methane  gas  emission  in  Pakistan.  Similarly,  energy 
consumption  Granger  cause  carbon  dioxide  emissions  (C02)  and 
energy  consumption  Granger  cause  combustible  renewables  and 
waste  (CRWASTE)  but  not  the  other  way  around.  However, 
agricultural  nitrous  oxide  emissions  (AGRNIT)  do  not  Granger 
cause  energy  consumption  and  energy  consumption  Granger 
cause  AGRNIT  is  accepted  via  both  route;  therefore,  we  may 
conclude  that  both  variables  are  causality  independent  in  nature. 
The  results  reflect  that  greenhouse  gases  are  closely  associated 
with  energy  consumption  and  GDP  per  unit  use  of  energy  in 


Pakistan.  In  reality,  the  agricultural  gases  are  closely  related  to 
energy  consumption  in  Pakistan. 


4.3.  Johansen  cointegration  test  between  all  greenhouse  gases  and 
energy  consumption  followed  by  variance  decomposition  analysis 

In  order  to  compare  the  contribution  extents  of  various  green¬ 
house  gas  emissions,  the  variance  decomposition  approach  is 
adopted  over  the  sample  period.  First,  we  take  the  energy 
consumption  as  the  dependent  variable,  while  greenhouse  gases 
coupled  with  GDP  per  unit  use  of  energy  as  independent  variables, 
and  conduct  the  Johansen  cointegration  test  among  these  vari¬ 
ables  over  a  period  of  36  years.  Table  6  shows  the  overall 
cointegration  results  of  greenhouse  gases  coupled  with  GDP  per 
unit  use  of  energy  and  energy  consumption. 

The  results  indicate  that  there  exists  statistically  significant 
cointegration  among  greenhouse  gases  and  energy  consumption 
in  Pakistan.  Next,  we  apply  the  variance  decomposition  approach 
based  on  the  vector  error  correction  model  (VECM)  to  explore  the 
influence  of  greenhouse  gases  on  energy  consumption,  and  com¬ 
pare  their  contribution  differences.  The  results  of  variance  decom¬ 
position  analysis  are  shown  in  Fig.  5. 

The  results  find  that,  among  all  greenhouse  gases,  Combustible 
renewables  and  waste  (CRWASTE)  in  the  form  of  metric  tons  of  oil 
equivalent  exerts  the  largest  influence,  whose  steady  contribution 
level  on  energy  use  to  36.35%;  subsequently  followed  the  influence 
of  agricultural  methane  emissions  (AGRMET)  in  the  form  of 
thousand  metric  tons  of  C02  equivalent;  GDP  per  unit  of  energy 
use  (GDPENRG)  in  the  form  of  constant  2005  PPP  $  per  kg  of  oil 
equivalent;  carbon  dioxide  emissions  (C02)  in  the  form  of  metric 
tons  per  capita  and  agricultural  nitrous  oxide  emissions  (AGRNIT) 


Table  4 

Results  of  Johansen  cointegration  tests. 


Cointegration  test 

Hypothesized  no.  of  CE  (s) 

Eigenvalue 

Trace  statistic 

5%  Critical  value 

Prob. 

Panel  A:  Series:  ENRG,  AGRMET,  GDPENRG 

None3 

0.456 

34.239 

29.797 

0.014 

At  most  1 

0.297 

12.901 

15.494 

0.118 

At  most  2 

0.015 

0.537 

3.841 

0.463 

Panel  B:  Series:  ENRG,  AGRNIT,  GDPENRG 

None3 

0.537 

36.795 

29.797 

0.006 

At  most  1 

0.250 

10.316 

15.494 

0.257 

At  most  2 

0.006 

0.226 

3.841 

0.634 

Panel  C:  Series:  ENRG,  C02,  GDPENRG 

None 

0.369 

26.017 

29.797 

0.128 

At  most  1 

0.240 

9.865 

15.494 

0.291 

At  most  2 

0.007 

0.254 

3.841 

0.613 

Panel  D:  Series:  ENRG,  CRWASTE,  GDPENRG 

None3 

0.481 

30.616 

29.797 

0.040 

At  most  1 

0.142 

7.619 

15.494 

0.506 

At  most  2 

0.061 

2.236 

3.841 

0.134 

Note:  dependent  variable  in  each  Johansen  cointegration  test  is  ENRG. 
a  Rejection  of  the  hypothesis  at  the  5%  level. 

Table  5 

Causality  test  results  among  greenhouse  gases  and  energy  consumption. 

Null  hypothesis 

Chi-square  statistic 

Prob. 

Decision 

ENRG  does  not  Grange  cause  the  changes  in  AGRMET 

8.223 

0.031 

Reject 

AGRMET  does  not  Granger  cause  the  changes  in  ENRG 

0.218 

0.896 

Accept 

ENRG  does  not  Grange  cause  the  changes  in  AGRNIT 

2.117 

0.346 

Accept 

AGRNIT  does  not  Granger  cause  the  changes  in  ENRG 

0.086 

0.957 

Accept 

ENRG  does  not  Grange  cause  the  changes  in  C02 

8.867 

0.019 

Reject 

C02  does  not  Granger  cause  the  changes  in  ENRG 

1.176 

0.554 

Accept 

ENRG  does  not  Grange  cause  the  changes  in  CRWASTE 

13.150 

0.001 

Reject 

CRWASTE  does  not  Granger  cause  the  changes  in  ENRG 

0.589 

0.744 

Accept 

Note:  The  modified  Granger  causality  test  approach  used  in  the  table  is  provided  by  Toda  and  Yamamoto’s  [38],  The  causality  tests  between  greenhouse  gases  and  energy 
consumption  are  based  on  the  significance  of  Chi-square  statistics  for  Wald  tests  of  VAR  models. 
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Table  6 

Johansen  cointegration  test  among  greenhouse  gases,  GDP  per  unit  use  of  energy  and  energy  consumption. 


Multivariate  cointegration  test 

Hypothesized  no.  of  CE  (s) 

Eigenvalue 

Trace  statistic 

5%  Critical  value 

Prob. 

Unrestricted  cointegration  rank  test  (Trace) 

None3 

0.695 

123.127 

95.753 

0.000 

At  most  l3 

0.640 

81.524 

69.818 

0.004 

At  most  2 

0.456 

45.711 

47.856 

0.078 

At  most  3 

0.292 

24.361 

29.797 

0.185 

At  most  4 

0.220 

12.258 

15.494 

0.144 

At  most  5 

0.096 

3.5612 

3.841 

0.059 

Unrestricted  cointegration  rank  test  (maximum  Eigenvalue) 

None3 

0.695 

41.602 

40.077 

0.033 

At  most  l3 

0.640 

35.813 

33.876 

0.029 

At  most  2 

0.456 

21.350 

27.584 

0.255 

At  most  3 

0.292 

12.102 

21.131 

0.537 

At  most  4 

0.220 

8.697 

14.264 

0.312 

At  most  5 

0.096 

3.561 

3.841 

0.059 

Note:  Trace  test  indicates  2  co-integrating  equations  at  the  0.05  level. 

a  Rejection  of  the  hypothesis  at  the  0.05  level.  Mackinnon  et  al.  [40]  p-values. 
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Fig.  5.  Variance  decomposition 
Source:  author’s  estimation. 

Note:  magnitude  of  the  coefficients  obtained  from  variance  decomposition  analysis. 

in  the  form  of  thousand  metric  tons  of  C02  equivalent,  with  steady 
contribution  level  of  5.21%;  4.71%;  1.89%  and  1.28%,  respectively. 
The  results  with  each  and  every  greenhouse  gas  emissions  with 
respect  to  their  influencing  magnitude  on  energy  consumption  are 
also  shown  in  Fig.  6  for  ready  reference. 


5.  Summary  and  conclusion 

The  objective  of  this  study  was  to  empirically  investigate  the 
effect  of  energy  consumption  on  greenhouse  gas  emissions  in 
Pakistan  over  the  period  of  1975-2011.  Overall,  the  analysis  shows 
that  there  is  a  long-term  relationship  among  energy  consumption, 
GDP  per  unit  of  energy  use  and  greenhouse  gas  emissions  in 
Pakistan.  Factors  such  as  the  agricultural  methane  emissions; 
agricultural  nitrous  oxide  emissions;  C02  emissions,  combustible 
renewables  and  waste  and  GDP  per  unit  of  energy  use  have  exerted 
the  contribution  of  5.21%;  1.28%;  1.89%;  36.35%  and  4.71%  changes 
in  energy  consumption  in  Pakistan.  As  far  as  the  results  of  Granger 
causality  test,  except  agricultural  nitrous  oxide  emissions;  energy 
consumption  Granger  causes  Greenhouse  gases  but  not  vice  versa. 
The  result  shows  unidirectional  causality  running  towards  energy 
consumption  to  Greenhouse  gas  emissions  in  Pakistan. 


Sustainable  environment  development  and  energy  efficiency 
are  essential  for  human  development.  The  national  challenges  of 
environmental  degradation  and  lack  of  access  to  clean,  affordable 
energy  services  also  have  a  global  impact,  as  climate  changes,  loss 
of  biodiversity  and  ozone  layer  depletion  cannot  be  addressed  by 
countries  independently  [29].  Chow  [30]  suggests  some  solutions 
related  to  the  problems  of  energy  and  environmental  degradation 
includes  (i)  reducing  the  use  of  energy  in  production  and  con¬ 
sumption,  (ii)  increasing  the  use  of  energy-saving  and  environ¬ 
mentally  friendly  methods  in  production  and  consumption  and 
(iii)  promoting  technological  innovations  that  will  reduce  the  use 
of  energy  per  unit  of  output  (reduce  energy  intensity  or  increase 
energy  efficiency)  or  reduce  pollution  per  unit  of  output. 

The  private  sector,  as  well  as  civil  society  at  large,  should  also 
be  encouraged  to  think  beyond  individual  or  corporate  interests 
towards  recognition  of  a  shared  responsibility  for  the  environ¬ 
ment.  Vigorous  resource  mobilization  to  combat  desertification 
would  stand  a  better  chance  of  succeeding  if  launched  on  the  basis 
of  empirically  verifiable  improvements  [31].  Education  of  the 
public  and  the  legislative,  political,  and  business  and  religious 
leaders  on  environmental  issues,  specifically  air  pollution  items 
and  greenhouse  gases  is  recommended.  Likewise,  communication 
and  coordination  between  environmental/health  and  safety/and 
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Variance  Decomposition  of  DLOG(ENRG) 
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Variance  Decomposition  of  DLOG(AGRNIT) 
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Variance  Decomposition  of  DLOG(CRWASTE) 
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Variance  Decomposition  of  DLOG(GDPENRG) 


DLOG(ENRG)  - DLOG(AGRMET) 

DLOG(AGRNIT)  - DL0G(C02) 

DLOG(CRWASTE)  - DLOG(GDPENRG) 


Fig.  6.  Variance  decomposition  analysis  with  other  greenhouse  gases  on  energy  consumption. 
Source:  author's  own  efforts  based  on  the  vector  error  correction  model  (VECM). 


other  agencies  is  highly  recommended  to  remove  barriers  to 
potential  mitigation  of  emissions  and  emission  effects  [32].  Due 
to  uncertainty  about  climate  change,  and  human  contributions 
thereto,  many  policymakers  call  for  "precautionary”  measures  to 
reduce  the  risk  of  global  warming.  Such  policies  are  characterized 
as  “insurance.”  Such  insurance  against  the  risks  of  climate  change 
can  be  achieved  by  either  lessening  the  likelihood  of  change  by 
reducing  atmospheric  concentrations  of  greenhouse  gases  through 
a  combination  of  emission  controls  and  carbon  sequestration 
strategies,  or  by  enacting  mitigation  measures  to  reduce  the 
possible  economic  and  ecological  impact  of  a  potential  climate 
change  [33].  Financial  incentives  can  be  helpful  to  kick-start  the 
market  for  new  energy  efficient  products  as  well  as  for  developing 
countries  where  funding  is  not  always  available  [34]. 

Climate  change  education  for  the  public  is  essential  to  inform 
rational  personal  choices.  Informal  science  education  can  play  a  critical 
role  in  reaching  broad  and  diverse  audiences  and  is  well  situated 
to  help  improve  public  climate  change  awareness,  understanding, 
and  informed  decision-making  [35].  The  challenge  of  managing  GHG 
emission  reductions  in  Pakistan  and  globally  is  significant,  creating 
both  immediate  and  long-term  research  needs.  Among  the  most 
important  near-term  research  needs  is  to  extend  the  economic 
analysis  conducted  here  to  additional  emission  reduction  strategies, 
and  to  improve  the  existing  modeling  capacity  for  advanced  energy 
systems.  Understating  the  gravity  of  the  situation,  we  need  to  put  such 
phenomena  under  the  straight  jacket  of  intuitional  control,  by  doing 
so;  we  can  lessen  the  loathed  impact  of  such  phenomena. 
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